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ABSTRACT 

We present new high spatial resolution HST/ ACS imaging of NGC 1140 and high 
spectral resolution VLT/UVES spectroscopy of its central star- forming region. The 
central region contains several clusters, the two brightest of which are clusters 1 and 
6 from Hunter, O'Connell & Gallagher, located within star-forming knots A and B, 
respectively. Nebular analysis indicates that the knots have an LMC-like metallicity of 
12 + logO/H = 8.29 ± 0.09. According to continuum subtracted Ha ACS imaging, 
cluster 1 dominates the nebular emission of the brighter knot A. Conversely, negligible 
nebular emission in knot B originates from cluster 6. Evolutionary synthesis modelling 
implies an age of 5 ± 1 Myr for cluster 1, from which a photometric mass of (1.1 ± 
0.3) X 10^ M0 is obtained. For this age and photometric mass, the modelling predicts 
the presence of ^ 5900 late O stars within cluster 1 . Wolf-Rayet features are observed 
in knot A, suggesting 550 late-type WN and 200 early-type WC stars. Therefore, 
N{WR)/N{0) ^ 0.1, assuming that all the WR stars are located within cluster 1. The 
velocity dispersions of the clusters were measured from constituent red supergiants as 
CT ~ 23 ± 1 kms^^ for cluster 1 and ct 26 ± 1 kms~^ for cluster 6. Combining a with 
half-light radii of 8 ± 2 pc and 6.0 ± 0.2 pc measured from the F625W ACS image 
implies virial masses of (10 ± 3) x 10^ Mq and (9.1 ± 0.8) x 10*^ Mq for clusters 1 
and 6, respectively. The most likely reason for the difference between the dynamical 
and photometric masses of cluster 1 is that the velocity dispersion of knot A is not 
due solely to cluster 1, as assumed, but has an additional component associated with 
cluster 2. 

Key words: galaxies: individual: NGC 1140 - galaxies: starburst - galaxies: star 
clusters - stars: Wolf-Rayet 



1 INTRODUCTION 

Violent bursts of star formation, which are characteristic of 
starburst galaxies, resemble the star-forming phase of young 
galaxies in the early Universe. Nearby starbursts provide lo- 
cal templates to which distant star-forming galaxies may 
be directly compared. Only a handful of starburst galax- 
ies are located within 10 Mpc, yet they produce around a 
quarter of the entire high-mass star population (|Heckmanl 
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Il998l ). Thus, starbursts are ideal sites in which to study mas- 
sive stars. The signatures of Wolf-Rayet (WR) stars, which 
are the highly evolved descendants of massive O stars, are 
apparent in a subset of starburst galaxies, known as WR 
galaxies. Since WR stars are exclusively associated with 
young (~ 5 Myr) stellar populations, WR galaxies repre- 
sent young or ongoing starbursts. A hallmark of all star- 
bursts seems to be the production of luminous, compact 
star clusters. The sizes, luminosities and masses of these 
Young Massive Clusters (YMCs) are consistent with the 
properties expected of young globular clusters (GCs). This 
led to the suggestion that YMCs represent globular clusters 
at an early phase of their evolution. Dynamical mass mea- 
surements can potentially be used to test the scenario with 
respect to the long-term survi vability of th ese YMCs (e.g. 
iHo fc Filippenkd Il996al lbl: see Ide Griis fc Parmentieil |2007| 
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for an overview). Understanding the role of YMCs as can- 
didate proto-GCs is vital for our understanding of galaxy 
formation and evolution, as well as large-scale star forma- 
tion. Clusters also have the advantage of being simple to 
model - they can be approximated as a coeval, simple stel- 
lar population with a single metallicity. This makes them 
ideal as probes of burst properties such as age, duration, 
chemical evolution and star-formation rate as well as con- 
straining the parameters of the stellar initial mass function 
(IMF). 

NGC 1140 is a low- metallicity WR galaxy at a dis- 
tance of ~ 20 and is a prime example of a nearby 
analogue of the star-forming galaxies identified at high 
redshifts. HI and optical observations indicate that the 
galaxy has undergone a merger within the past 1 Gyr 
and it is thought that this event is responsible for the 
plethora of YMCs that the galaxy hosts. Star clusters 
have been imaged with both the Wide Field/Planetary 
Camera (WF/PC) and the Wide Field Planetary Cam- 
era 2 (WFPC2) on the Hubble Space T elesc ope (HST) by 
iHunter. O'Connell fc Gallagher fl994a1l and lde Griis et al.l 
l|2004h . They identified eight young luminous clusters in the 
centre of the galajcy, and their studies indicated that the 
clusters have masses of up to a few x 10® M© and ages that 
lie in the range of a few to a few tens of Myr. 

This paper considers the properties of the two bright- 
est clusters within NGC 1140, clusters 1 and 6, and the two 
star-forming knots in which these clusters are contained. It 
is structured as follows. Details of our VLT and HST obser- 
vations and data reduction are given in Section 2, along with 
photometry of the clusters. In Sections 3 and 4, we discuss 
the stellar and nebular properties of the two knots apparent 
in the VLT spectra, and in Section 5 we present the results 
of evolutionary synthesis modelling of cluster 1. The mas- 
sive star population of cluster 1 and of both knots A and B, 
and the star- formation rate of NGC 1140, are considered in 
Section 6. The dynamical masses of the clusters are deter- 
mined in Section 7 and we discuss our findings in Section 8. 
Finally, we summarise our results in Section 9. 



2 OBSERVATIONS AND DATA REDUCTION 

We obtained high-resolution spectroscopy of the central re- 
gion of NGC 1140 with the VLT/Uhraviolet and Visual 
Echelle Spectrograph (UVES), in addition to high spatial 
resolution Hubble Space Telescope (HST) / Advanced Cam- 
era for Surveys (ACS) imaging of NGC 1140. Fig. [T] shows 
the region of the galaxy observe d with UVES and lab els 
clusters 1 — 7, as designated by iHunter et al. fl994a')l . It 
shows the location of the 1x11 arcsec^ UVES slit, which was 
aligned north-south over the two knots in the central region 
of NGC 1140, hereafter called knots A and B. Knot A, which 
contains the clusters 1 and 2, is the brightest region in the 
optical and near infrared. Knot B lies ~ 3 arcsec south of 
knot A, and hosts clusters 5, 6 and 7. 



Table 1. F625W ACS aperture STMAG photometry of the bright- 
est central clusters of NGC 1140 and the regions designated knots 
A and B. The F300W a nd F814W cl uster magnitudes are taken 
from WFPC2 imaging of lde GriiseTal . (2004). 



Cluster 




mpsoow 


'TIF625W 


mF814W 






(mag) 


(mag) 


(mag) 


Cluster 


1 


16.50 ± 0.05 


17.9 ± 0.1 


18.51 ± 0.05 


Cluster 2 


16.37 ± 0.05 


18.8 ± 0.1 


19.11 ± 0.08 


Cluster 


6 


18.47 ± 0.09 


18.4 ± 0.1 


18.73 ± 0.03 


Cluster 


7 


21.71 ± 1.51 


19.7 ± 0.1 


19.86 ± 0.04 


Knot A 




15.5 ± 0.1 


17.0 ± 0.1 


17.8 ± 0.1 


Knot B 




17.6 ± 0.1 


17.5 ± 0.1 


18.4 ± 0.1 



2.1 Imaging of NGC 1140 

NGC 1140 was observed on 21 August 2003 with ACS/WFC 
aboard HST as part of GO programme 9892 (PI Jansen). 
Two exposures were taken - one of 950s duration with the 
narrow-band F658N filter and one of 100s duration using 
the broad-band F625W filter. The data were reduced using 
the on-the-fly reduction pipeline. 

At the distance of NGC 1140, the F658N fiher does 
not include the [N II] 6583 line, and, therefore, is essentially 
an Ha plus continuum filter. We have corrected for the un- 
derlying continuum by subtracting the scaled F625W image 
such that a background galaxy is removed. This continuum 
subtracted F658N is included in Fig. [1] The high resolu- 
tion of the ACS data also shows that clusters 3 and 5 from 
iHunter et al. (1994a)l are not single clusters. 

Aperture photometry was carried out on the drizzled 
F625W ACS image in Starlink's GAIA. Since the crowded 
nature of the field makes it difficult to position apertures 
to include all of the light from the desired cluster without 
contamination from its neighbours, several circular aper- 
tures with a range of radii and positions were considered 
for each region. Each measurement was sky subtracted us- 
ing a mean sky value obtained from several small apertures 
of sky, and aperture c orrected using the values contained in 
ISirianni et al. (2005)1 . As these aperture corrections are valid 
for point sources rather than for spatially extended clusters, 
our photometry may be an underestimate. However, since 
the light profile of the clusters is difficult to ascertain (see 
Section [7. ip . no correction for this effect is made. 

Table [l] presents our F625W photometry. The associ- 
ated uncertainty represents the range of results produced. 
Table [T] also contains F300W and F814W photometry. 
The F300W and F814W photometry of knots A and B 
were determined using archival WFPC2 images (GO 8645, 
PI Windhorst) using the same method as described for 
the F625W images. The F300W and F814W c luster pho- 
tometry was taken from Ide Griis et al. (|2004 ). The pre- 
COSTAR WF/PC ph otometry of lde Griis et al.1 (|2004l l and 
IHunter et al. (1994a)l was not included due to its lower spa- 
tial resolution. 



^ Based on the heliocentric velocity of the galaxy, cor- 
rected for the Virgocentric flow, and assuming Hq = 
70 km s~^Mpc~^ ; adopted from the HyperLeda database at 
|http://leda .univ-lyonl.fr/l 



2.2 Spectroscopy of NGC 1140 

The central region of NGC 1140 was observed with UVES 
on the VLT Kueyen Telescope (UT2) in Chile on 19-20 
September 2003. UVES is a cross-dispersed echelle spec- 
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Figure 1. HST / ACS images of NGC 1140, at difTcrent cut levels (1 arcscc ~ 100 pc), marked as being either the F625W image or the 
continuum subtracted F658N image. The left-hand image shows the whole galaxy. The white box indicates the region that is magnified 
in the images in the right-hand panel. The lighter diagonal band across the image is the inter-CCD gap between the two detectors of 
ACS. Superimposed on each right-hand image is the 1 X 11 arcsec^ VLT/U VES slit. The region s extracted as knots A and B and as sky 
are indicated by the solid and dashed lines, respectively. Clusters 1 — 7 (iHunter et al. 1994al ) are also labelled. The top right-hand 
image shows that clusters 3 and 5 likely comprise several smaller clusters. 



trograph, with two arms - a blue arm comprising a single 
EEV CCD and a red arm comprising a mosaic of an EEV 
and a MIT-LL CCD. Using the dichroic beam splitter, red 
and blue data were taken simultaneously: cross-disperser #4 
with 312 I mm~^ was centred on 840 nm and cross-disperser 
#2 with 660Zmm~^ was centred on 437 nm. Thus, data 
from the regions ~ 3740 - 4985A and ~ 6475 - 10090A were 
obtained. The spectral region ~ 8180 - 8420A was not ob- 
served since it lay on the gap between the two red CCDs. 
The data were 2x2 binned, so that the red data have a 
pixel scale of 0.364 arcsec pixel"^ and the blue data have a 
scale of 0.512 arcsec pixel"^. 

The galaxy was observed with UVES for a total time 
of 14700s on the first night and 16200s on the second night. 
The seeing averaged ~ 0.8 arcsec on the first night, varying 
by ~ 0.1 arcsec and was somewhat poorer on the second 
night. The two knots were more clearly resolved in the data 
obtained on 19 September and so only the data of NGC 1140 
from this night were considered in the analysis. 

A telluric star, HIP 11918, was also observed in the 
same setup. Flux standards were not part of the original 
observing programme. However, we obtained data of two 
flux standards, EG274 and EG21. EG274 was observed for 
another programme on the first half of the first night for 
600s using the blue setup above, plus a red setup centred 
on 8600 A. EG21 was observed on the second night in 4 ex- 
posures of 2s in only the blue setup. ThAr arcs were taken 
for wavelength calibration. The resolution of the data, as 
measured from the ThAr arc lines, is ~ 7kms~^. 

The data were reduced using the software package iraf. 
The data were bias corrected with a median bias frame, and 
divided by a normalised flat field. The two knots were ex- 
tracted, each with an aperture of ~ 1.3 arcsec, and back- 
ground subtracted using a polynomial fit to the sky regions 
indicated in Fig. [T] After extraction, median knot spectra 
were created and wavelength calibrated, before the orders 
were merged. The standard stars were similarly extracted. 



This method over-subtracts the sky in the region of neb- 
ular lines. Therefore, the NGC 1140 spectra were also ex- 
tracted without sky subtraction. The sky level was then 
manually subtracted with a polynomial fit so that the con- 
tinuum matched that of the sky-subtracted spectra away 
from the nebular lines. 

The blue NGC 1140 data were fiux calibrated with 
EG274, and adjusted to match the F625W photometry of 
knot A and B contained in Table [1] Due to the lack of a 
suitable red fiux standard, it was not possible to flux cali- 
brate the red data. The calibrated blue spectrum of knot A 
is shown in Fig. [2j with the principal nebular lines indicated. 

2.3 Spectroscopy of red supergiants 

Comparison red supergiant (RSG) spectra are required to 
determine a cluster's velocity dispersion based upon its RSG 
features, which can then be used to calculate its dynamical 
mass. Therefore, six Small Magellanic Cloud (SMC) red su- 
pergiants, taken from thCiMassey &__Qlscn ( 2 003. ) catalogue, 
were observed with the setup described above. The catalogue 
ID number of these stars, and their spectral types were 20133 
(MO I), 64448 (K2-7 I), 71566 (K7 1), 50840 (Ml-2 I), 66694 
(K5 I) and 71507 (K3-5 I). 



3 STELLAR PROPERTIES 

Red supergiant and OB stellar features are visible in the 
UVES spectra of both knots. Wolf-Rayet (WR) features are 
also present in the spectrum of knot A. 

3.1 OB and WR star features 

Knots A and B both show stellar absorption in the hydrogen 
Balmer series, arising from the presence of early-type stars, 
underlying strong nebular emission. The equivalent widths 
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Figure 2. Blue spectrum of knot A with the principal nebular 
lines and the Wolf-Rayet blue bump marked, corrected to the 
stellar rest- frame. Broad stellar absorption is apparent underlying 
the narrow nebular Balmer emission. 



of the underlying stellar absorption are Wa(H/9) ~ 4A, 
Wa(H7) ~ 3A and Wa(H5) ~ SA. 

WR stars are the bare cores of O stars, in their final 
evolutionary stages, characterised by strong, broad emis- 
sion lines. Their strong winds reveal first the CNO-burning 
products, observable as nitrogen-rich (WN) stars and sub- 
sequently the He burning products, producing carbon- 
rich (WC) stars. A WR galaxy, or a cluster containing 
WR stars, can be identified by the presence of the blue 
bump around A4686A or the ye llow bump at A5808A 
(|Schaerer. Contini fc Pindao 1999a ). The yellow bump is 
produced solely by WC stars, while the blue bump can con- 
tain features of WN stars, such as stellar N III 4640 and He II 
4686, and/or WC stars, such as C HI 4650 and C IV 4658, 
with narrow nebular features superimposed, the strongest of 
which is [Fe HI] 4658. 

The blue WR bump emission feature in the spectrum 
of knot A has an equivalent width of Wa ~ 12. 9A, with 
He II 4686 contributing around ~ 4.6A. The presence of 
WR features in the spect rum of knot A constrains the age 
of the knot to 4 ± 1 Myr l|Crowtheijl2"007l '). No WR emission 
features are visible in the spectrum of knot B. The absence of 
WR stars implies that knot B is either younger than around 
3 Myr, or older than around 5 Myr. However, the presence 
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Figure 3. The normalised Ca II triplet lines of knot A (top) and 
knot B (bottom) after correction for a recession velocity of Vr = 
1475kms-l. The A8542 and the A862 lines are offset by 0.5 and 
1.0, respectively. Nebular Paschen emission lines are indicated. 



of RSG features in the spectrum (see Section [32} makes it 
unlikely that the knot is younger than 3 Myr, for the case 
of an instantaneous burst. 



3.2 RSG features 

The Ca II triplet absorption lines (AA8498, 8542, 8662 A) are 
clearly visible in the red spectra of both knots A and B. 
These features, along with other weak metal lines present in 
the data, such as Fe I and Mg I, arise from the RSGs in the 
cluster. The presence of RSGs indicates a cluster age greater 
than around 5 Myr. 

The Ca II triplet lines of both knots were fitted with 
Gaussian profiles with the elf (emission-line fitting) rou- 
tine in Starlink's dipso package. The line centres, widths 
and intensity of the fits were allowed to vary freely, and the 
measured line centres used to determine the recession ve- 
locity of the galaxy. A good fit could not be obtained to 
the A8542 and A8662 lines of knot A, due to the presence 
of Paschen nebular lines within these profiles. Therefore, the 
recession velocities determined from the remaining four lines 
were averaged to obtain Vr ~ 1475 ±4kms~^. Our result is 



Cluster and nebular properties of NGC II40 5 



in good agreement with other recent values of recessional ve- 
locity, such as Vr = 1498 ± 33kms ~^ from optical measure- 
ments (|de Vaucouleurs et al.l 199lh. Vr- = 1501 ± lkms~^ 
from 21-cm HI observations ( Havnes et al.|[l"998[) and Vr = 
1480 ± gkms^'^ from [Fe II] measurements~ (|de Griis et alj 
l2004h . 

After velocity correction, the equivalent widths of the 
Ca II triplet lines were measured with elf, fixing the central 
wavelengths of the lines. This yields values of Wa(8498) = 
0.77 ± O.O7A, Wa(8542) = 1.05 ± 0.2lA and Wa(8662) = 
0.91 ± O.ISA for knot A and Wa(8498) = 0.92 ± 0.09A, 
Wa(8542) = 2.2 ± O.SA and Wa(8662) = 1.45 ± 0.15A for 
knot B. The velocity corrected Ca II triplet lines for knots 
A and B are presented in Fig. |3l 




-100 -50 50 
Velocity (km s"^ ) 



4 NEBULAR PROPERTIES 

In this section we derive information on the gas dynamics of 
the knots from the profiles of the nebular lines. Extinctions, 
electron densities, temperatures and elemental abundances 
are determined from the relative fluxes of the nebular lines. 



4.1 Line profiles and dynamics 

Fig. m shows that the nebular H/3 profiles of knot A and 
knot B are very different. Knot A is dominated by one com- 
ponent at approximately v ~ — 15kms~^, but also has a 
weaker underlying broad component. Knot B more clearly 
comprises two components: a weaker component blueshifted 
at V ^ -36kms"^ with FWHM ~ 58kms~^ and a brighter 
redshifted component at « 28kms~^ with FWHM ~ 
47kms~^. These nebular profiles are representative of all 
the strong nebular lines of the knots. Examination of the 2D 
image indicates that there is diffuse emission throughout the 
central region of the galaxy that probably corresponds to the 
bluer, weaker component of knot B and to the weak, broad 
component of knot A. The stronger components of the knots 
arise from more discrete emission from the clusters within 
these knots, as are visible in Fig.[T] Alternatively, the broad 
underlying component of the knot A profiles may be due to 
hot cluster winds impacting on the surrounding interstellar 
clou ds, as has been se en in other young star-forming regions 
(seelSidoU et al.ll2006l for a review and IWestmoguette et al.l 
I2OO7 for a detailed discussion) . Turbulent broadening of the 
order of 30 km s~^ was implied from the FWHMs of the main 
components of the strong nebular lines of knots A and B. 



4.2 Line fluxes and extinction 

Nebular line fluxes were measured by fitting the nebular 
emission with Gaussian profiles using the elf routine, al- 
lowing line centres, widths and intensity to vary freely. The 
strong nebular emission of both knots A and B was fitted 
with two Gaussians (Section |4. If) . 

The internal extinctions of the knots, E{B — V)int, 
were determined by comparing the observed fiux ratios of 
F(H7)/F(H^) and F(H^)/F(H/3) to the intrinsic intensity 
ratios predicted by assuming Case B recombination theory 
for electron densities of 10 ^ cm"'' and a temperature of 10* K 
IIHummer fc Storevlll987l ). As the difi'use emission contained 



Figure 4. The H/3 profiles for knots A and B, after correction to 
Vr = 1475 km s~^. 



Table 2. Comparison between predicted intrinsic flux ratios, 
Ia/I(H/3) and observed flux ratios, Fa/F{H/3). For the fluxes 
quoted here, both components of the knot A profile were summed, 
while only the strong redshifted component of knot B was consid- 
ered. The observed ratios were first dereddened with a standard 
Galactic extinction law. Internal extinctions, E{B — V )int, were 
then d erived from predicted intensity ratios assuming a lHowarthl 
lll983l) LMC extinction law. 





Line 


Fa/F(H/3) 


Ia/I(H/3) 


E{B - V)int 






(Observed) 


(Predicted) 


(mag) 


Knot A 


H7 


0.427 ± 0.021 


0.468 


0.16 ± 0.10 




H<5 


0.225 ± 0.011 


0.259 


0.16 ± 0.07 


Knot B 


H7 


0.427 ± 0.021 


0.468 


0.19 ± 0.10 




H<5 


0.212 ± 0.011 


0.259 


0.27 ± 0.07 



within knot B is expected to have a very low extinction com- 
pared to the discrete emission from the clusters, the nebu- 
lar extinction was calculated based only on the stronger, 
redder component of knot B. Since the division between 
components is less well defined for knot A, both compo- 
nents were summed. The observed flux ratios were first cor- 
recte d for foreground extinction of E(B — y)fore = 0.038 
mag (|Schlegel. Finkbeiner fc Davislll998l ) using a standard 
Galactic extinction law before determining the internal ex- 
tinctions of the knots using a L arge Magellanic Cloud (LMC) 
extinction law (iHowarthll 19831 ). The ratios measured and in- 
ternal extinction calculated are presented in Table [2] Ratios 
of H7 and H(5 to H/3 were considered since Ha was not flux 
cafibrated. Mean values of E{B — V)int = 0.16 ± 0.06 mag 
for knot A and E{B - V)int = 0.23 ± 0.08 mag for knot B 
were subsequently adopted. 

Table |3] presents observed and extinction corrected line 
fluxes, normalised to H/3 = 100. Here, both components were 
summed for both knots, as the total emission within the 
knot was of interest. Lines with total fluxes of < 0.2% of H/3 
were not included. We concentrate on deriving the electron 
density, temperature and nebular abundances for knot A be- 
cause of the added complications introduced by the nebular 
structure present in Knot B. 
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Table 3. Observed (Fa) and intrinsic nebular line fluxes rel- 
ative to H/3 = 100 for knot A and knot B. The fluxes quoted 
here are the summation of both components of the nebular lines, 
for both knots A and B. The intrinsic values are corrected for a 
foreground extinction of E{B — V)f^^^ = 0.038 mag with a stan- 
dard Galactic extinction law and then for internal extinctions of 
E{B - V)int = 0.16 ± 0.06 mag for k not A and E(B - V)int = 
0.23 ± 0.08 mag for knot B using the iHowarthI lll983h LMC ex- 
tinction law. The flux of H/3 is in units of xlO~^* ergs~^ cm~^. 



Table 4. Equivalent widths and knot ages implied from the 
Starburst99 model. The e quivalent width of Ca II is define d as 
Wa{A8542) -I- Wa(A8662) Jpiaz. Terlevich fc Terlevichlll989l) . 



Knot A 



Knot B 



Line 




Fx 


h 


Fx 


Ix 


A3726 


[O II] 


65±3 


82±7 


81±4 


112±12 


A3729 


[O II] 


89±4 


lllitlO 


110±6 


151±17 


A3869 


[Ne III] 


21. Oil. 


25.5±2.0 


20.0±1.0 


26.1±2.5 


A3967 


[Ne III] 


5.7±0.6 


6.8±0.8 


4.8±0.5 


6.1±0.8 


A4026 


He I 


0.89±0.18 


1.05±0.22 


0.89±0.18 


1.11±0.23 


A4069 


[S II] 


0.95±0.19 


1.10±0.23 


2.27±0.23 


2.8±0.3 


A4102 


H<5 


22.6±1.1 


26. Oil. 7 


23.8±1.2 


28.9±2.2 


A4341 


H7 


42.7±2.1 


46.8±2.7 


43.1±2.2 


49±3 


A4363 


[O HI] 


1.5±0.15 


1.63±0.17 


1.03±0.10 


1.16±0.12 


A4471 


He I 


3.1±0.3 


3.3±0.3 


3.4±0.3 


3.7±0.4 


A4658 


[Fe HI] 


1.08±0.11 


1.12±0.11 


2.22±0.22 


2.33±0.23 


A4686 


He II 


0.59±0.12 


0.61±0.12 


0.88±0.18 


0.91±0.18 


A4702 


[Fe HI] 


0.52±0.10 


0.53±0.11 






A4861 


H/3 


100 


100 


100 


100 


A4959 


[O HI] 


98±5 


96±5 


96±5 


94±5 


H/3 




5.4±0.5 


10±3 


2.70±0.27 


7±3 



4.3 Electron densities and temperatures 



The electron density, A^e 



Te, of 



knot A were determined in the five-level atom calculator 
TEMDEN within IRAF using the diagnostic line ratios of [O II] 
I(A3726)/I(A3729) and [O III] I(A4363)/I(A4959). This gives 
values of A'^e = 60 ± 50 cm~^ (consistent with the low density 
limit) and Te = 9700 ± 300 K. 



4.4 Abundances 

Knot A abundances were calculated from standard [O II] 
and [O III] diagnostics, plus the values of A^e and Te deter- 
mined in Section|331 These yield values of iV(0+)/A''(H+) = 
(8.2 ± 2.0) X 10"^ and iV(0^+)/iV(H+) = (1.10 ± 
0.20) X 10~^. These imply an abundance for knot A of 
12 -I- logO/H = 8.29 ± 0.09. We adopt this value for 
NGC 1140. It lies between the abundance of the Small Mag- 
ellanic Cloud (SMC) of 12 + log O/H = 8.13 and that o f 
the LMC of 12 + log O/H = 8.37 l|RusseU fc Dopitalll990D . 
It agrees well with other dire ct abundance measurem ents 
of NGC 11 40 of 8.18 ± 0.06 Jlzotov fc ThuanI 120041 ') and 
8.26 ± 0.07 (|Nagao. Maiolino fc Marconill2006ir^ ecalculated 
from the measurements of Izotov fc Thuanll2004l ') . 

As the determined abundance is most similar to the 
LMC abundance, we adopt an LMC-like metallicity with 
Z — 0.008 for spectral modelling purposes. 





Knot A 


Knot B 


Line 


Wx{A) 


Age (Myr) 


Wx{k) 


Age (Myr) 


Ha 


289±29 


5.5±0.5 


142±14 


6.0±0.5 


H/3 


57±6 


5.5±0.5 


62±6 


5.0±0.5 


Ca II 


1.96±0.28 


5.5±0.5 


3.7±0.3 


>6 



5 PHOTOMETRICALLY DETERMINED 
CLUSTER PROPERTIES 

Having established that knot A has an LMC-like metallic- 
ity and an age of around 5 Myr, we computed v5.1 Star- 
burst99 evolutionary synthe sis models (|Leitherer et al.|[l999l : 
IVazauez fc Leithererl I2OO5I ) for an instantaneous burst of 
star formation, with a total stellar mass of 10® Mq and a 
metallicity of Z = 0.008 for age s betwee n 1 an d 10 Myr in 
0.5 Myr intervals. We adopted a lKroupal (|2002l ) IMF (with 
slope a = 2.3 for 0.5 ^ M/Mq ^ 100 and a = 1.3 for 
0.1 M/Mp) < 0.5) and Padova stellar evolutionary tracks 
(|Fagotto et al.|[T994l ). which include careful consideration of 
the RSG phase. 



5.1 Age 

The Starburst99 model predicts the equivalent widths for 
Ha, H/3 and the Ca II triplet. The equivalent widths mea- 
sured for each knot and the ages implied from the Star- 
burst99 model are presented in Table ID 

Fig. [1] and Table [1] help to put these equivalent widths 
in context by showing the regions in which the continuum 
light and the nebular emission originate. In knot A, both the 
stellar continuum and nebular emission are dominated by 
cluster 1, indicating that the Ha and H/3 equivalent widths 
are meaningful and can be applied to cluster 1. On the other 
hand, the stellar continuum of knot B is dominated by clus- 
ter 6, while the nebular emission is due solely to cluster 5 and 
a second source of emission ~ 0.4 arcsec immediately to the 
southwest of cluster 5. Therefore, these latter two sources 
must be younger than the 5 — 6 Myr quoted in Table ID 
Indeed, the lack of Wolf-Rayet features in knot B indicates 
that they are younger than ~ 3 Myr. Fig. [1] shows no neb- 
ular emission from cluster 6 or 7, so these clusters must be 
older than ~ 10 Myr, unless any gas associated with these 
clusters has been removed by multiple supernovae at an even 
earlier stage. The equivalent widths of the Ca II triplet for 
both knots A and B are meaningful and can be applied to 
cluster 1 and 6, respectively, since these clusters dominate 
the continuum in the region. However, only a lower limit to 
the age of cluster 6 can be determined from the equivalent 
width of the Ca II triplet, because it does not vary strongly 
with age after a few Myr. Furthermore, the contribution 
of RSGs to cluster ev olution is not well understood (e.g. 
iMassev fc Olsenll2003l ). and as such RSGs do not provide a 
reliable cluster age indicator. 

After considering both the age indication of 4 ± 1 Myr 
from the presence of WR features in the spectrum of knot A 
and the realistic age estimates from Table |4l we adopt an 
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Wavelength (A) 

Figure 5. Fit of tlie reddened 5 Myr old model spectrum (Star- 
burstgg) to ACS and WFPC2 photometry for cluster 1. The 
model spectrum was first redd ened by the foregrou nd extinction 
of E{B - V)iore = 0.038 mag jSchlegel et aljliggsl ) using a stan- 
dard Galactic r eddening law and then further reddened with a 
(iHowarthlllisi ) LMC extinction law to best fit the plotted pho- 
tometry. 



5.3 Photometric mass 

The F625W magnitude for cluster 1 was converted into 
an apparent T/-band m agnitude using the transformation 
in ISirianni et al. (2005)1 . assuming the Starburst99 model 
V — R colour for a 5 Myr old burst and adopting the clus- 
ter extinction determined in Section 15.21 A distance of 20 
Mpc to NGC 1140 then implies a V-band luminosity of 
Lv ~ (7.3 ± 1.8) X 10^ Lq for cluster 1. Comparing this 
luminosity to the Starburst99 prediction yields a photomet- 
ric mass of (1.1 ± 0.3) X 10^ M0. This is based on a Kroupa 
IMF. 



6 MASSIVE STAR POPULATION AND 

STAR-FORMATION RATE OF NGC 1140 

In this section we determine the O star content of the knots 
from the H/3 luminosity of the UVES spectroscopy. We es- 
timate the WR content of knot A from the blue bump seen 
in its spectrum and measure the star formation rate (SFR) 
of NGC 1140 from the continuum subtracted F658N ACS 
imaging. 



age of 5 ± 1 Myr for knot A. Since no age can be estimated 
for cluster 6, it is not considered further in this section. 

We find an age difi'erential betwe en knot A and knot B, 
in a greement with t he res ults of both iHunter et al. fl994all 
and Ide Griis et al.1 (l2004h . Based on their V — I colours, 
iHunter et al. (1994a)l estimated the ages of the clusters in 
the region of knot A as ~ 3 Myr and interpreted the redder 
colours of the knot B clusters as being due to an older age of 
these clusters, as opposed to a higher extinction. They esti- 
mated the ages of these clusters as ~ 15 Myr, since this was 
the ag e determined by lO'ConneU. GaUagher III fc Hunteii 
l| 19941 ) for the YMCs in NGC 1569 and NGC 1705 
wi th the same co l ours. The near infrared CIRPASS data 
of Ide Griis et al.l l|2004l ) showed strong [Fe II] emission 
throughout the galaxy, while strong Br(12-4) and Br(14- 
4) was predominantly confined to knot A. Since the diffuse 
[Fe II] emission is likely associated with supernova remnants, 
while the Br ackett nebular emissio n lines are associated with 
HII regions, Ide Griis et al.l (|2004l ) argued that the ratio of 
these lines is a good age indicator. This would suggest that 
knot B is several Myr older than knot A. 



5.2 Extinction 

The model spectrum produced by the Starburst99 model for 
a 5 Myr old instantaneous burst was reddened to reproduce 
the ACS and WFPC2 photometry of cluster 1 (Table [T|. 
The best fit, presented in Fig. (5] yields a value of inter- 
nal cluster reddening of E{B — V)int ~ 0.24 ± 0.04 mag 
for cluster 1. The uncertainty quoted here solely considers 
the photometric uncertainties. The F300W and F625W pho- 
tometry was more heavily weighted than the F814W pho- 
tometry, since these are more sen sitive to redden ing. Our 
value agrees well with the result of Ide Griis et al.1 (2004) of 
E{B - V)int = 0.14 - 0.26 mag for cluster 2, determined 
from spectral energy distribution fits to their WF/PC and 
WFPC2 photometry. It also agess with the nebular extinc- 
tion determined from the Balmer lines, given in Tabled 



6.1 O star content 

An estimate of the number of O stars within an individ- 
ual knot can be obtained from the H/3 luminosity, L(H/?). 
Assu ming Case B recombination theory ([Hummer fc StorevI 
Il987l ) and that a 'normal' stellar population is the only 
source of ionising photons within the knot, the number of 
equivalent 07V stars, A'^(07V), contained within the knot 
is given by: 



N{07Y) 



LjUp) X 2.1 X 10^^ 



(see e.g. IVaccal Il994l ). Here, Q()^^ is the observed total 
Lyman continuum luminosity of the knot and Q'o^^ is 
the Lyman continuum fiux of an individual 07V star, 
distance of 20 Mpc and using the flux values given in 
Table O we obtain L(H^) = (5.0 ± 1.8) x 10^^ erg s"^ 
for knot A and L(H/3) = (3.1 ± 1.9) x lO^^ergs"^ for 
knot B. Taking Qp '^^ = 10'^ ^-^ photons s"\ as suggested by 
i Hadfield fc Crowtherl l|2006l ) for LMC metallicity O stars, 
implies that knot A contains 1300 ± 500 equivalent 07V 
stars and that knot B contains 800 ± 500 equivalent 07V 
stars. The Starburst99 modelling predicts an H/3 luminos- 
ity for cluster 1 of (4.5 ± 1.1) x lO^^ergs"^, showing that 
cluster 1 dominates knot A in terms of its H/3 luminosity. 
This implies the presence of 1200 ±300 equivalent 07V stars 
within cluster 1. 

This number of equivalent 07V stars can be converted 
into the total number of O stars, N{0), using the time- 
dependent parameter, ri{t): 



iV(0) 



iV(07V) 



For an age of 5 Myr and a metallicity of Z — 0.008, 
ISchaerer fc Vaccal j 19981 ) give rj = 0.2, implying that clus- 
ter 1 contains 5900 ± 1400 O stars. As knot A is dominated 
by cluster 1, applying rj = 0.2 provides a good approxima- 
tion to the number of O stars within knot A of 6600 ± 2400 
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relative contributions to the H/3 flux of the central region of 
2:1 (taken from Table|3]). Using the equation: 



SFR(MQyr- 



7.9 X 10 
1 



X 10"^^ X 



X L(Ha)(ergs"^) 

^ (photons s^^) 

iKennicuttl liggsl ). we estimate the SFR of the 
galaxy of 0.7 ± O.3M0yr~^, in good agreement 
with the SFR of 0.8MQyr"\ obtained from 
ground-based continuum subtracted H a imaging by 
iHunter. van Woerden fc Gallagher ('1994b)l . after adjust- 
ment to our adopted distance. Our Ha luminosity gives 
gObs ^ (-g g _j_ 2.6) X lO''^ photons s-i or 9000 ± 3000 07V 
equivalents in NGC 1140. 



Figure 6. The A4640 Wolf-Rayet blue bump in tlie spectrum of 
knot A. The solid line shows the best empirical fit to the spectrum, 
of 550 WNL stars and 200 WCE stars. The dashed line shows 
only that from 550 WNL stars and the dotted line shows only the 
contribution from 200 WCE stars. Nebular [Fo III] A4658 emission 
is marked. 

O stars. For knot B, which is dominated by the young clus- 
ter 5, the total number of O stars is comparable to the num- 
ber of equivalent 07V, since ?7 ~ 1 for ages of < 3 Myr 
l|Schaerer fc Vaccalll998l ). 

6.2 Wolf-Rayet star content 

The number of WR stars within knot A can be obtained 
by considering the luminosity of the blue WR bump. The 
similarity in the strengths of N HI 4640 and He H 4686 
in knot A suggests the presence of late-type WN (WNL) 
stars. The luminosity of A4686 is 6.8 x 10^®ergs"\ after 
correcting for an extinction of E{B — V)int ~ 0.24 mag and 
a distance of 20 Mpc. Comparing this to the average lu- 
minosity of LMC WN7-9 st ars of (7.2 ± 6.7) x lO^'^ergs"^ 
ICrowther fc HadfieldllioO^ . indicates that knot A contains 
950 WNL stars, assuming no contribution to the blue bump 
from WC stars. However, iGuseva. Isotov fc ThuanI (|2000l ) 
identify early-type WC (WCE) stars in their long-slit spec- 
trum of the entire central region of NGC 1140, by the pres- 
ence of C IV 5808 emission. Therefore, we empirically match 
the blue bump using various multiples of an LMC WN7-9 
star, plus an LMC WC4 star, finding a best fit of approx- 
imately 550 ± 100 WNL stars and 200 ± 50 WCE stars. 
This is shown in Fig. (6] We assume that all of the WR stars 
within knot A are contained in cluster 1, since it is the most 
massive young cluster within this region. 

6.3 Star formation rate of NGC 1140 

We measure an observed Ha fiux of F(Ha) = (1.2 ± 0.2) x 
lO"^'^ ergs~^ cm~^ for the whole galaxy, based upon aper- 
ture photometry of the continuum subtracted F658N image 
(see Section [2TT|) . Assuming a distance of 20 Mpc and a mean 
extinction of E{B — V)int ~ 0.19 mag, we find an Ha lu- 
minosity of NGC 1140 of (9 ± 3) X 10*° ergs"^ This mean 
extinction was determined as an average of the nebular ex- 
tinctions of knots A and B (Section 14. 2|) weighted by the 



7 DYNAMICAL MASSES OF THE CLUSTERS 

The virial motions of the individual stars within a cluster, 
as measured by the line-of-sight velocity dispersion, give a 
dynamical estimate of the cluster mass. This also requires a 
knowledge of the half-light radius of the cluster, which is the 
two-dimensional radius within which half of the projected 
luminosity of the cluster is contained. 

7.1 Half-light radius 

The half-light radii of clusters 1 and 6 were determined 
from the driz zled A CS/ HST F625W images using the pro- 
gram ISHAPE (|Larsen 1999). The routine convolves the point- 
spread function (PSF) of the camera with the desired light 
profile. ISHAPE computes a minimisation between the 
cluster image and the model profile over a set fitting radius, 
iterating the FWHM, the ratio of the minor to major axes 
and the orientation of the fit. In order to obtain the best 
fit, the user specifies the 'clean radius', which is the largest 
radius that contains no contamination by neighbouring clus- 
ters. 

Cluster 6 is more isolated than cluster 1, and is uncon- 
taminated by its neighbours up to a clean radius of four 
pixels (1 pixel — 0.05 arcsec). Cluster 1, however, is only 
isolated up to two pixels. The half-light radii were com- 
puted by ISHAPE for a range of profiles, using a P SF that 
was computed in Tiny Tim (|Krist fc Hooklll997h . for fit- 
ting radii between two and twenty pixels with a clean ra- 
dius of two pixels for cluster 1 and four pixels for cluster 6. 
The profiles considered wore a Gaussian profile, EFF pro- 
files (Els on. Fall fc Freeman. ,1987 ) of index x, which take 
the form: 



and iKingj (|1962| ) profiles of index c of the form: 

f(z)=\ (vi+^" ^'^'^ 

[ c. 

There was a sudden jump in the values of half-light radii 
produced for fitting radii of eight pixels and for nine pixels. 
Therefore, the means of the half-light radii computed by 
ISHAPE for each profile for fitting radii between two and eight 
pixels were considered. These, along with the corresponding 
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Table 5. Mean half-light radii computed by ISHAPE for a range 
of model profiles for fitting radii between two and eight pixels, 
and the corresponding standard deviations (s.d.). Half-light radii 
were adopted on the basis of the best standard deviation. 



Cluster 1 Cluster 6 



Model 


Index 


fhi (pc) 


s.d. 


rill (pc) 


s.d. 


Gaussian 




6.1 


0.8 


5.6 


0.4 


EFF 


15 


9.5 


0.5 


8.4 


0.5 


EFF 


25 


6.9 


0.8 


6.12 


0.12 


King 


5 


6.6 


0.8 


5.79 


0.12 


King 


15 


8.2 


0.5 


7.0 


0.6 


King 


30 


10.4 


0.4 


8.8 


1.2 


King 


100 


17.3 


0.6 


14.6 


2.6 


Adopted 




8 ± 2 




6.0 ± 0.2 



standard deviation, are contained in Table [5] As a good pro- 
file fit will not vary much with fitting radius, the profile with 
the lowest standard deviation was adopted. For cluster 1 this 
was the King 30 profile. However, the standard deviation of 
this is quite large and similar to that produced for other pro- 
files. Discarding the outlying values produced by the King 
100 profile, the mean and standard deviation of the 42 re- 
sults produced by the other 6 profiles gives rhi = 8 ± 2 pc 
for cluster 1. 

For cluster 6, both the EFF 25 and the King 5 profile 
show an equally high level of consistency and a mean of the 
results of both of these profiles was adopted to give ru = 
6.0 ± 0.2 pc for cluster 6. 

7.2 Velocity dispersion 

The line-of-sight velocity dispersion, a, of a cluster can be 
determined by comparing the broadening in the lines of the 
cluster spectrum with respect to those in a red supergiant 
template spectrum. There are two main methods that quan- 
tify the comparison, and these are discussed below. 

The spectrum is compared to red supergiant templates 
because the lines of these stars are broadened by only 
a few kms~^, by m acro-turbulence in their atmospheres 
l|Grav fc Toned 1 19861 ) . It is not appropriate to use ear- 
lier type supergiants and main sequence stars, as effects 
such as rotational broadening, macro-turbulence and micro- 
turbulence broaden the lines of these stars by amounts 
comparable to the anticipated cluster velocity dispersions. 
Therefore, only spectral regions redwards of ~ 5000A, 
which are dominated by light arising from cool super- 
giants of spectral type F -M, should be considered (see e.g. 
IHo fc Filippenkq Il996ah . Suitable spectral regions should 
show visible similarity between the cluster and RSG spectra 
and should also not contain any telluric lines, which provide 
an artificial match between the spectra. 

The first method relies on minimising a reduced 
between the normalised cluster spectrum and normalised, 
broadened template spectra. The broadening is achieved by 
convolving the normalised template spectrum with a Gaus- 
sian of a equal to the desired velocity broadening. The 
broadened template is multiplied by an optimum factor, so 
that it produces the lowest possible reduced-x^ between the 
broadened template spectrum and the cluster spectrum. The 



Table 6. Values of line-of sight velocity dispersion over the re- 
gion 8485 — 8845A for each red su pergiant template obse rved. 
The templates are identified by their iMassev fc OlsenI (|2003| ) cat- 
alogue number. It is marked whether the cross-correlation (Xcor) 
or minimisation (x^) technique was used. The Paschcn emis- 
sion was masked out for the reduced x^ minimisation technique. 
The results obtained for the cross-correlation technique are an 
upper limit on the cluster velocity dispersion. 



RSG Template Velocity dispersion (kms ^) 



Cat No. 


Spectral 


Clusterl 


cluster 6 




Type 


X^ 


Xcor x^ 


64448 


K2-7 I 


23 


26 26 


30840 


K3-5 I 


23 


29 27 


66694 


K5 I 


23 


30 26 


71566 


K7 I 


22 


29 26 


20133 


MO I 


24 


31 27 


50840 


Ml-2 I 


22 


32 26 


Mean 




24 


30 26 



template is broadened by a range of suitable values and a 
reduced-x^ is similarly computed for each value of broaden- 
ing. The lowest of these x^ values determines what amount 
of broadening produces the best match with the cluster for 
a given template, and so indicates the velocity dispersion of 
the cluster. This is repeated for a range of template spec- 
tral types. This method relies on a good match between the 
relative line strengths of the cluster and template and can 
thus be sensitive to the spectral type of the template star. 

The second method utilis es the cross-correlation tech- 
nique of lTonrv fc David (|l979D . It requires that the spectra 
being considered are normalised, and continuum subtracted, 
to give a flat continuum at zero. The cluster is then cross- 
correlated with a red supergiant template over suitable spec- 
tral regions, and the FWHM of the resulting cross correla- 
tion function (CCF) is measured. The template spectrum is 
broadened by a range of velocities. Each broadened template 
is cross-correlated with the original, unbroadened template 
and the FWHM of each CCF is measured. In this man- 
ner, the near-linear relationship between broadening and the 
FWHM of the CCF can be empirically calibrated to an ab- 
solute scale. This calibration is applied to the FWHM of the 
original CCF, produced by cross-correlating the cluster spec- 
trum with the template spectrum, to determine the velocity 
dispersion of the cluster. This is repeated for each template 
star. The conversion factor differs with template spectrum. 
While this method is less sensitive to spectral type matching, 
it suffers from complications associated with the subjectiv- 
ity of fitting CCFs. These include factors such as selecting 
the background level and fitting non-Gaussian CCFs, and 
are especially important when the CCF is weak. 

Since clusters 1 and 6 are the brightest cluster members 
in knots A and B, we assume that the velocity dispersions of 
these clusters dominate the broadening of the RSG features 
apparent in the knots. 

Unfortunately, the youth of cluster 1 means that the 
RSG features in the spectrum are very weak, and many of 
the lines visible in the template spectra are absent in the 
cluster spectrum. The low signal-to-noise ratio of knot B 
causes a similar problem. The strongest RSG features, the 
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Ca II triplet absorption lines, are clearly visible in the cluster 
spectra. However, these lines are saturated in the template 
spectra. As the core of a saturated profile is narrower than 
would have been produced in a Gaussian profile, and it is the 
core that produces the CCF signal, cross-correlation of these 
regions t ends to overestimate the cluster velocity dispersion 
(see e.g. lWalcher et al.ll2005h . However, cross-correlation of 
all other regions that both contain RSG features and lack 
telluric features produce very noisy, non-Gaussian CCFs, 
which cannot be robustly fitted. Therefore, only the region 
8485 - 8845A was considered for cross-correlation. Cross- 
correlation of even this region of knot A with the tem- 
plates did not produce CCFs that could be confidently fitted, 
likely due to the contamination of the Ca II triplet by the 
strong Paschen lines. Therefore, no results were produced 
for knot A using the cross-correlation technique. 

The reduced-x^ minimisation could also only be com- 
puted for the region 8485 — 8845 A, due to the lack of strong 
RSG absorption lines except around the Ca II triplet. The 
Paschen emission features were masked out of both knot 
spectra before the reduced-x'^ minimisation was computed. 
The results of the cross-correlation and the reduced-x^ min- 
imisation are listed in Table |6] 

The results show consistency over all the spectral types, 
with the cross-correlation of cluster 6 with the templates 
producing systematically higher velocity dispersions than 
the reduced-x^ technique, as expected. There should not 
be any systematic uncertainties in the velocity dispersion 
calculated by reduced- x'^ minimisation, des pite the use of 
the Ca II triplet lines. Mengel et al.l ()2002l 'l found no dis- 
parity between the velocity dispersion results computed by 
X^ minimisation for the strongest component of the Ca II 
triplet and other individual absorption features for clus- 
ters in NGC 4038/4039. The mean and standard devia- 
tion of the velocity dispersions calculated from the reduced- 



Table 7. Summary of cluster properties of clusters 1 and 6. 



o- = 24 ± Ikms" 
cluster 6. 



for cluster 1 and cr = 26 ± 1 km s for 



T.3 Virial mass 

The virial equation relates the virial mass of a cluster, Mdyn, 
to the line-of-sight velocity dispersion, a, and the half-light 
radius, rhi, of the cluster by the equation: 



dyn 



G 



This equation assumes that a cluster is gravitationally 
bound, spherically symmetrical and virialised, that the ve- 
locity dispersion of the cluster is isotropic and that all of the 
stars containe d in the cluster are single stars of equal mass. 
ISpitzeil (|l987l ) showed that r) = 9.75 for globular clusters 
with a wide range of light profiles. However, mass segrega- 
tion and the presence of binary stars within the cluster can 
cause a large var iation in the value of r; (Boily et al. 2005; 
iFleck et al. 20061 : Kouwenhoven & de Grijs, in prep.). For 
the young age and high mass of cluster 1, however, these 
effects are minimal. Table [7] summarises the cluster proper- 
ties. 



Property 



cluster 1 cluster i 



E{B - y){clustcr) (mag) 


0.24±0.04 




Age (Myr) 


5±1 




rill (pc) 


8±2 


6.0±0.2 


(T (km s~^) 


24±1 


26±1 


Mdyn(lOSM0) 


10±3 


9.1±0.8 




7.3±1.8 




Lv/MiyMLv/M)Q) 


7.0±2.6 





8 DISCUSSION 

In this section we discuss the disparity found between the 
dynamical and photometric mass estimates. We also con- 
sider the reliability of stellar evolution models by comparing 
them with observations. 



8.1 How reliable are photometric mass estimates 
versus dynamical mass estimates? 

The dynamical mass determined for cluster 1 is several times 
greater than the photometric mass found. To assess this dis- 
parity. Fig. [7] compares the light-to-dynamical mass ratio 
of the cluster to a model light-to-mass r atio at the known 
age of the cluster. The model assumes iMarastonI (|2005l ) 
SSPs, a Kroupa IMF and solar metallicity. Data for sev- 
eral other clusters whose velocity dispersions were measured 
from UVES observations are also included. It is clear from 
this figure that many other young (< 20 Myr) clusters also 
have light-to-dynamical mass ratios well below model pre- 
dictions, while the older (> 20 Myr) clusters have ratios that 
agree well with the canonica l value. 

As discussed by iBastian et al. I l|2006t ) and 
ICoodwin fc BastianI (|2006l ). the most likely explana- 
tion for the discrepancy between model and observations 
seen in the sample of young clusters is a lack of virial 
equilibrium caused by violent relaxation after the formation 
of the cluster and from the expulsion of gas. This causes 
the measured dynamical mass to overestimate the true 
mass. Since the older clusters, which are expected to be 
in virial equilibrium after surviving gas expulsion, lie on 
or near the model line, the photometric mass is likely to 
be a good representation of the true cluster mass. This 
supports the assumption of a standard IMF. In this case, 
the disagreement between the dynamical mass and the 
photometric mass of a cluster can be used to assess to what 
degree the cluster is out of virial equilibrium. This can 
be parameterised by the effective star formation efficiency 
(eSFE) of the cluster, e; at the onset of gas expulsion, a clus- 
ter with an eSFE — e has a velocity dispersion t hat is 
too large to be in virial equilibrium (Goodwin &i BastianI 
2006). The model light-to- mass ratio assuming a range of 
effective star formation efficiencies are included in Fig. [71 
taking the onset of gas expuls ion as 2 Myr. These tracks 
differ from those presented by iGoodwin fc BastianI (|2006l ) 
as they are not smoothed. 

An alternative explanation for clusters lying below the 
canonical line in Fig. [7] is that these clusters have non- 
standard IMFs. In this situation, the dynamical mass repre- 
sents the true cluster mass, and the incorrect assumption of a 
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Kroupa IMF 



eSFE 
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Figure 7. Diagram of light to dynamical mass ratio against 
age for cluster 1 in NGC 1140 (filled triangle; this work), 
and for clusters 1 2 and 16 in NGC 4038/4039 (filled circles; 
iMengel et al] 120021). clusters 502 and 805 in N GC 5236 (open 
trian gles; iLarse n fc Richtler||2 004h. clusters W3 ijMaraston et al] 
\2004) and W30 llBastian et al.l2006h in NGC 7252 (fille d squares), 
cluster G114 in NGC 1316 (cross: iBastian et aIll20C)6h arid clus- 
ters #23 and #34 in ESQ 338-IG04 (open circles: lOstlin et alj 
l2006h . Also plotted are the model evolution of light-to-mass ra- 
tio predicted by a Kroupa IMF with Maraston {2005} SSPs for 
a solar metallicity burst and the effect that non-100% effective 
star-format ion efficiencies (eSFEs) h ave on the model light-to- 
mass ratio ijCoodwin fc Bastianll2006h . 



standard IMF causes the photometric mass to differ from the 
dynamical mass. Since the discrepant clusters lie below the 
model line, they would have an over-abundance of low-mass 
stars with respect to the Kroupa IMF. Such bottom-heavy 
clusters are more likely to remain bound than clusters with 
a standard Kroupa IMF and should be apparent in the sam- 
ple of old clusters. As they are not, this explanation seems 
unlikely for the majority of these young clusters. However, 
a non-standard IMF can not be ruled out in the case of 
NGC 1140-1. 

There are various potential explanations for the differ- 
ence between light-to-dynamical mass ratio and model light- 
to- mass ratio in the case of NGC 1140-1. These are explored 
below: 

(i) Assuming that the difference bewteen dynamical and 
photmetric mass is due solely to the cluster being out of 
virial equilibrium, it can be see n from Fig. [7] that cluster 1 
has an eSFE of around 10-20%. [Goodwin fc BastianI (120061 ) 
predict that a cluster with an eSFE of 10% will lose 80% 
of its mass in ~ 20 Myr and not remain bound. As such, 
the cluster would disperse over a relatively short timeframe 
and, therefore, not evolve into a second generation globular 
cluster. 

(ii) An underestimate of the cluster extinction would have 
caused an underestimation of the photometric mass. How- 



ever, in order to increase the luminosity of the cluster by 
the factor of ~ 9 required to match the photometric mass 
with the dynamical mass, the extinction would need to have 
been underestimated by around E{B — V)int ~ 0.9 mag. 
Since the extinction calculated is in good agreement with 
previous measurements, and is already larger than the neb- 
ular extinction derived, this is unlikely to be the case. 

(iii) Determining the half-light radii of the clusters was 
heavily dependent on the hght profile adopted and several of 
the profiles were in good agreement with observation. How- 
ever, adopting rhi ~ 6.0 pc, the lowest potential half-light 
radius computed by ishape for cluster 1, would still pro- 
duce a virial mass that was a factor of seven larger than the 
photometric mass. 

(iv) The crowded nature of the region may have caused 
the cluster velocity dispersion to be overestimated. The 
broadening of the RSG features seen in the knot spectra 
may not have been due solely to the virial motions of the 
brightest cluster within the knot, as is assumed, but may 
have had contributions from the motions of other clusters. 
This situation can be modelled using the template spectra, 
and a reduced minimisation can be computed for this 
model knot. A template spectrum is broadened by 8kms~^ 
to represent the broadening of cluster 1, as implied by its 
photometric mass and half-light radius. Added to this is the 
cluster 2 contribution, in the fiux ratio of 1/1.7, as implied 
by the F814W photometry in Table [T] Taking the velocity 
dispersion of cluster 2 to be equal to that of cluster 1, the 
velocity dispersion found for knot A in Table [6] can be re- 
produced by introducing a relative velocity shift between 
the two clusters of ~ 20kms~^. It is not possible to assess 
how likely such a shift is. Cluster 2 appears to be older than 
cluster 1 because of the lack of nebular emission (Fig. [1]). 
This then implies that cluster 2 could have a larger pho- 
tometric mass than cluster 1, despite its fainter luminosity, 
and thus have a larger velocity dispersion than adopted here. 
This would reduce the velocity shift required to produce a 
combined velocity dispersion equivalent to that observed for 
knot A. However, neither the photometric mass nor the shift 
are known and the importance of the contribution of cluster 
2 to the velocity dispersion of knot A cannot be assessed. 

It is not possible to quantify the degree to which these four 
factors contribute to the difference between the light-to- 
dynamical mass ratio measured and the model in Fig. [T] 
However, the contribution from the neighbours of cluster 1 
is likely to be the most dominant factor for the large dynam- 
ical mass measured. As it is not possible to quantify this, it 
is also not possible to assess how long the cluster is likely to 



8.2 Young clusters as testing grounds of stellar 
evolution models 

The nebular properties of the two knots of NGC 1140, and 
the massive star content of clusters 1 and 6 are presented in 
Table[Sl There have been several other studies of young mas- 
sive stellar populations of nearby starbursts, including the 
ultrav iolet (UV) survey of IChandar. Leitherer fc Tremontil 
(j2004l) and optical studies of Ivacca fc Cont il (|1992| ) and 
I Schaerer. Contini fc Kunth (1999a)l . However, the results 
of these surveys are not readily comparable to our results. 
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Table 8. Nebular knot properties and massive star content of NGC 1140 (this work), NGC 3125 l|Hadfield 
dSidoli et al.ll2006l) . 



Crowthejr2006l ') and Tol 89 



Galaxy 


NGC 1140 


NGC 3125 


Tol 89 (NGC 5398) 


Knot 


A B 


A B 


A B 



E{B - V) (mag) 


0.16±0.06 


0.23±0.08 


0.16 


0.13 


0.07 


0.29 


Afe (cm-3) 


60±50 




140 


140 


90 ±40 


150 ± 40 


T. (K) 


9700±300 




10500 


9800 


10000 ± 300 


9800 ± 300 


12 + log O/H 


8.29±0.09 




8.32 


8.35 


Q Q7+0.08 
°-^'-0.09 


a qO+0.06 


N{or 


6600±2400 


800±500 


4000 


3200 


685 


2780 


Age (Myr) 


5±1 


> 10 


~ 4 


~ 4 


3.5-4 


< 2.5 



Cluster 


1 


6 


Al 


A2 


Bl + B2 


Al 


Bl 


Mphot(lO«M0) 


1.1±0.3 




0.20 


0.22 


0.16 


0.1 - 0.2 


~ 0.03 


Ar(0)t 


5900±1400 




550 


600 


450 


310 - 680 


~ 100 


Af(WN) 


550±100 




105 


~ 55 


40 


80 




Af(WC) 


200±50 




20 




20 







N(WR)/N(0) 


0.1 




0.2 


0.1 


0.1 


0.2 





* From the luminosity of the nebular H/3 emission of the knot. 

t From Starburst99, for the age and photometric mass of the cluster. 



due to the different techniques adopted, especially in deter- 
mining the WR populations. 

IChandar (|2004l ) derived the O and WR content of 
the central starburst regions for a large sample of galaxies 
based upon UV HST spectroscopy and Starburst99 mod- 
els. Fits to far-UV spectral morphologies and continuum 
slopes provided representative ages and extinctions from 
which O star populations were obtained, although individ- 
ual clusters were not distinguished. The number of WR stars 
was inferred from a calibration of He II 1640 line lumi- 
nosities. Individual clusters wer e also not considered i n the 
ground-based optical s tudies of IVacca fc Coiitil (Il992l ) and 
ISchaerer et al. fl999a)l . Standard nebular techniques were 
apphed to derive extinctions, representative ages and O 
star numbers. WR populatio ns resulted from calibrations of 
Galactic and LMC WR stars (jVaccall 1994 1 Schaerer fc Vaccal 
1 19981 ). which differ from our calibration with solely LMC 
templates. 

Studies of individual massive clusters using common 
tec hniques, including HST ima ging, have been carried out 
bv lHadfield fc Crowtheij (|2006l ) and ISidoli et"aLl (|2006l ) for 
two galaxies of LMC-like metallicity, the blue compact dwarf 
galaxy NGC 3125 and the giant HII region (GHR) Tol 89 
located in the barred spiral galaxy NGC 5398. A comparison 
of individual knots, as derived from nebular properties, and 
cluster masses and O star numbers, from Starburst99 mod- 
elling, is made in Table [S] WR populations foll ow from opti- 
cal calibrations of ICrowther fc Hadfieldl (|2006l ) in all cases. 

NGC 1140-1 is a significantly more massive counter- 
part to NGC 3125-Al and Tol 89-Al, with a similar, high 
ratio of WR to O type stars, A''(WR)/iV(0) ~ 0.1 - 0.2 and 
iV(WC)/Af(WN) ^ 0.4. These empirical stellar populations 
can be compared to predictions from evolutionary synthe- 
sis models (e.g. Starburst99), for which reduced WR popu- 
lations of N{WR)/N{0) ~ 0.02 - 0.1 and a substantially 
different subtype distribution of iV(WC)/iV(WN) ~ 40 
are predicted. These are based on standard Geneva evo- 
lution models llMeviiet et all Il994l ) for ages of 4 — 5 Myr 
(see also ISidoh et al.ll2006l ). Allowance for rotational mix- 



ing in evolutionary models does not resolve these discrepan- 
cies, i.e. N( WR)/N(0) < . 1 and 7V(WC)/iV(WN) ~ 10 
according to IVazguez et al.l (|2007l ) for a 5 Myr old pop- 
ulatio n at Z = 0.008. Padova models (|Fagotto et al.l 
Il994h also predict a low ratio of WR to O stars, with 
A''(WR) /A''(0) ~ 0.04 - 0.07. These models, however, pre- 
dict significantly more realistic ratios of WC to WN stars of 
iV(WC)/iV(WN) ~ 0.2 - 2 between 4-5 Myr. 



9 SUMMARY 

We present new high spectral resolution VLT/UVES spec- 
troscopy and HST / ACS imaging of the central region of 
NGC 1140. It is apparent from the ACS imaging that the 
central region contains several clusters, although this is only 
resolved into two star-forming knots by the UVES spec- 
troscopy: knot A, which contains clusters 1 and 2, and 
knot B, which contains clusters 5, 6 and 7. 

Nebular analysis of knot A yields an LMC-like metal- 
licity of 12 + log O/H = 8.29 ± 0.09. Starburst99 mod- 
elling indicates an age of 5 ± 1 Myr and a photometric 
mass of (1.1 ± 0.3) X 10*^ M0 for cluster 1. Virial masses 
of (10 ± 3) X lO'' M© and (9.1 ± 0.8) x 10^ M© were deter- 
mined for clusters 1 and 6, respectively, using the half-light 
radii determined from the F625W ACS image and the ve- 
locity dispersions determined by computing a reduced-x^ 
minimisation between the cluster spectrum and RSG tem- 
plate spectra over a spectral region containing the Ca II 
triplet. We interpret the difference between the dynamical 
and photometric mass of cluster 1 as due to the crowded 
nature of knot A: the velocity dispersion measured may not 
relate only to cluster 1, as assumed, but likely contains a 
component that arises from cluster 2, with a different sys- 
temic velocity to cluster 1. 

We find 6600 and 800 O stars within knots A and 
B, respectively, from the H/3 luminosities of the knots. 
Our Starburst99 model predicts 5900 O stars within clus- 
ter 1. This implies that > 90% of the O stars within 
knot A are contained in cluster 1. Empricial fitting of the 
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blue bump of cluster 1 indicates that this cluster contains 
around 550 WN stars and 200 WC stars, giving ratios of 
iV(WR)/iV(0) = 0.1, if all of the WR stars lie within clus- 
ter 1, and iV(WC)/iV(WN) = 0.4. The observed ratio of WR 
stars to O stars is significantly larger than predicted by cur- 
rent evolutionary models. The observed ratio of WC to WN 
stars is reproduced more successfully using Padova instanta- 
neous burst models than Geneva models, even allowing for 
rotational mixing. 
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